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In a weak measurement with post-selection, a measurement value, called the weak value, can be 
amplified beyond the eigenvalues of the observable. However, there are some controversies whether 
the weak value amplification is practically useful or not in increasing sensitivity of the measurement 
in which fundamental quantum noise dominates. In this paper, we investigate the sensitivity limit 
of an optical interferometer by properly taking account quantum shot noise and radiation pressure 
noise. To do so, we formulate the weak value amplification in the Heisenberg picture, which enables 
us to intuitively understand what happens when the measurement outcome is post-selected and the 
weak value is amplified. As a result, we found that the sensitivity limit is given by the standard 
quantum limit that is the same as in a standard interferometry. We also discuss a way to circumvent 
the standard quantum limit. 


I. INTRODUCTION 


The idea of weak-value amplification (WVA) was orig¬ 
inally introduced by Aharonov, Albert, and Vaidman 
(AAV) in 1988 [I| (see 0 for a review). When a system 
is weakly measured by a measuring device, the measure¬ 
ment results can be much larger than the eigenvalues of 
the observable by appropriately selecting initial and final 
states of the system. This theoretical prediction has been 
demonstrated in various pioneering experiments, e.g. the 
rotation of photon polarization [3,13: quantum box prob¬ 
lem [3], the arrival time of a single photon [3], the spin 
Hall effect of light Q , optical beam deflection [3,0 , and 
optical phase [^. For further recent theoretical and ex¬ 
perimental developments on the WVA, see review papers 

The interesting nature of the WVA results from the 
definition of the weak value: 




( 1 ) 


where A is an observable associated with the system to 
be measured, |'0i) and \'4>f) are the initial and final states 
of the system. The weak value is interpreted as the ob¬ 
servable evaluated at intermediate times between the pre- 
and post-selections. From this definition, we see that if 
the pre- and post-selected states are nearly orthogonal, 
it seems that the weak value becomes arbitrarily large. 
In fact, the experiments mentioned above have shown 
that the WVA significantly improves signal-to-noise ratio 
(SNR) in the situation where technical noise (e.g. align¬ 
ment noise) dominates and outperforms a standard in¬ 
terferometry. The theoretical study also confirms this 
advantage of the WVA, e.g. [iliil: though it has been 
pointed out that the applicability of the amplification 
strongly depends on the property of technical noise [l3| . 
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On the other hand, when the weak value becomes large, 
we must take into account the nonlinear effects of the 
von-Neumann measurement [TtI - I^ . As the result, the 
amplification factor has a maximum value and vanishes 
when the pre- and post-selected states are exactly or¬ 
thogonal. Even in the weak measurement regime, the 
nearly orthogonal post-selection severely reduces output 
statistics, which consequently compensates the improve¬ 
ment of sensitivity due to the amplification of the signal. 
On the other hand, some authors [U, [23] have claimed 
that the sensitivity is infinitely improved by optimizing 
the wave functions of a system and a probe. However, if 
quantum shot noise in the detection process is taken into 
consideration, their conclusion would change. Therefore, 
the practical usefulness of the WVA in increasing sensi¬ 
tivity has been often controversial. 

In an optical interferometry, it was pointed out in 
that if photon shot noise dominates, there is no 
advantage to use the WVA for improving the fundamen¬ 
tal limit of parameter estimation precision. From an in¬ 
formational approach with quantum Fisher information 
[ 234 ^ . it is also obtained the same conclusion that the 
WVA is suboptimal in enhancing the estimation precision 
and does not perform better than the standard statisti¬ 
cal techniques (for the controversy on this conclusion, 
see also [23 - 1^ ). Indeed, these conclusions have been 
demonstrated in an optical experiment [s^- However, 
in an optical interferometry, not only shot noise but also 
radiation pressure noise contribute to fundamental quan¬ 
tum noise and leads to a kind of sensitivity limit, the 
so-called standard quantum limit (SQL) [^ |^. It is 
not clear how radiation pressure noise affects the sensi¬ 
tivity limit of the measurement and alters the SQL from 
that in the standard interferometry when the WVA is 
implemented. These questions are worth investigating 
because there are many systems, in which the radiation 
pressure and the SQL pl^ an important role, such as 
optomechanical systems [33l l34| and gravitational-wave 
detection [sslls^j. 

Conventionally, the weak measurement with post- 
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selection and the WVA are formulated in the Schrddinger 
picture. However, this is somewhat less intuitive, com¬ 
pared with the Heisenberg formulation, and often leads 
to confusion. Also it is more difficult to deal with 
moving mirrors and the radiation pressure noise in the 
Schrddinger picture. Instead, here we formulate the 
WVA in the Heisenberg picture and deal with the quan¬ 
tum noises in a fully quantum-mechanical way, including 
the radiation pressure noise. Then we derive the SQL 
and discuss a method to circumvent beyond the SQL. 

This paper is organized as follows. In Sec.|lTl we briefly 
review position measurements in an optical interferome¬ 
ter and the WVA in the Schrddinger picture. In Sec. IHIl 
we formulate the WVA in the Heisenberg picture, ap¬ 
propriately taking into account quantum shot noise and 
radiation pressure noise. The SQL is derived in Sec. 113 
and the method to overcome the SQL is presented in 
Sec. 13 Finally, Sec. ED is devoted to a summary. In this 
paper, we use the unit c = I. 


II. WEAK VALUE AMPLIFICAITON OF A 
SIGNAL IN THE SCHRODINGER PICTURE 

We briefly review the Schrddinger formulation of WVA, 
which is the original formulation by Aharonov et al. [l| 
and conventionally used so far. In this paper, we concen¬ 
trate on a small displacement measurement of mirrors 
in an optical interferometer, particularly focusing on a 
Michelson interferometer shown in Fig. [TJ This does not 
loose generality because a Michelson interferometer is ge¬ 
ometrically equivalent to a Mach-Zehnder interferometer 
and a Sagnac interferometer and gives the same sensitiv¬ 
ity limit of a position measurement. 



Figure 1. Michelson interferometer. 


We start with a single photon case for an illustrative 
purpose. (We will generalize the result to a macroscopic 
beam later.) A photon is injected from the left side, 
which we call the bright port because we will extend later 


to the case of a laser light pulse. The photon that en¬ 
ters the interferometer takes two paths with the same 
unperturbed distance L after being divided by the beam 
splitter. One of the photons (beams) in the y-arm is 
phase-shifted by 6 (pre-selection), reflected at the mirror 
with sensing small displacements i of the mirror (weak 
interaction), and recombined at the beam splitter with 
another photon that sensed small mirror displacement 
in the x-arm (post-selection). Here we are interested 
in small differential displacement of the mirrors since we 
focus on the detection of external disturbances on the 
mirrors. Finally, the photon is detected by the photode¬ 
tector at the output (dark) port. 

In this optical setup, we regard the beam’s which-path 
information (y-arm or x-arm) as the system to be mea¬ 
sured via a weak measurement. The initial pre-selected 
state of the system is denoted by 

The initial phase offset 6 is symmetrized merely for sim¬ 
plicity of calculation. An observable of an measuring 
device, or a pointer variable, is photon’s momentum (fre¬ 
quency), which measures the phase shift induced by the 
mirror displacements in the Michelson interferometer. 
The initial state of the pointer is 

1$) = J dp$(p)|p) . 

Since we measure the small displacement of the mir¬ 
rors at the asymmetric output port of this optical con¬ 
figuration, the observable is 2£A, where the operator 
A= |y)(y| — I®)(a;I carries information about which arm 
a photon passes and has the eigenvalues ±1, depending 
on the photon path. 

The Hamiltonian of the interaction, which is switched 
on at time to, is written as 

H = g5{t — to) A® p . (2) 

Here we defined g = —21. This interaction Hamiltonian is 
interpreted as a generalization of von Neumann interac¬ 
tion like £®p to the Michelson-type interferometer. After 
the interaction given in Eq. m and the post-selection by 
the final state of the system, \^|;f) = {\y) — \x))/^/2, the fi¬ 
nal state of the device can be ex actly evaluated including 
nonlinear terms in the coupling [23|, 

|d>') = 

= J dp‘^{p)\p){'tpf\tjji){cosgp-lA^singp) , (3) 

The expectation value of the n-th power of p for the 
final state of the photon in Eq. © is given by 

, ^ (p") + (|A^P - I)(p" sin^ gp) +lmA.uj{p'^ sm2gp) 
^ 1 -k (I 12 - 1) (sin^ gp) + ImA^, (sin 2gp) 

( 4 ) 
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The bracket (• • •) and (• • •) denote averaging over the 
initial and final state of the measuring device, respec¬ 
tively. In our case, the weak value defined in Eq. (HD 
is 

0 

= -i cot - . (5) 

The series of power in Eq. ([4]) contains all information 
about the system. If we measure the shift of the pointer 
variable (p) , the variance Var[p] = (p^) — ((p) is 
regarded as frequency noise, because it results from a 
photon spectral distribution. However, in an optical ex¬ 
periment, shot noise coming from the fluctuation of the 
photon number also contributes. As shown in [^ . taking 
into account the photon number fluctuations, it turned 
out to be that the shot noise is given by (p^) and is 
always larger than the frequency noise. Therefore, it is 
essential to consider the shot noise when we discuss sen¬ 
sitivity limit in an optical interferometer. 

Hereafter we denote p by fiw as we are considering an 
optical experiment. Given the initial momentum distri¬ 
bution of a photon is non-zero-mean Gaussian (For mul¬ 
tiple photons, a pulsed laser whose central frequency is 
mode-locked to wg), 


$(u;) = 


27rCT2 ) 


1/4 


exp 


(w - UJoY 


4ct2 


and substituting Eq. ([5]) for Eq. (|4]) , we obtain the follow¬ 
ing expressions for the first and second powers of w — wg 


0 : 


g{uj-UJo) = 


g^iico-LOor) 


s sin(0 -|- 4>) 

1 — cos(0 -I- (()) 

2s cos(0 -I- (/)) 


1 + 


where 


s = = 8ay , 


1 — cos{d + (j))_ 

(j) = —2gujQ = 4ujg£ . 


( 6 ) 

(7) 

( 8 ) 


The parameter s characterizes the measurement 
strength, since large g means strong coupling of the in¬ 
teraction and large <7^^ means the narrow distribution of 
photons in the time domain. Then the weak measure¬ 
ment is defined by the limit, s —>■ 0, which corresponds to 
two physical situations: (i) the measurement coupling, in 
this case mirror displacement, is small, (ii) probe’s wave 
function is narrow in momentum (frequency) space or 
broad in spatial (time) domain. 

At the limit of weak measurement, if only linear terms 
in s are kept, the above equations, and © , are re¬ 
duced to 

g{bj — Wg) ~ S COt(0 -I- 4)) , 

g^{{uj-uj,f)' (9) 


^ The sign of 0 is different from that in [2^ . This is just the matter 
of a different sign for mirror displacement. 


From these expressions, one finds that the frequency shift 
is proportional to cot(0 -|- (p) and could be amplified for 
small (j) by taking small 6. Note, however, that these ex¬ 
pressions are valid only in the linear regime of s and that 
sensitivity is not infinitely amplified if nonlinear orders 
are included. 

The above results can be easily extended to multiple- 
photons. From Eq. o, the probability distribution for a 
single photon at the output is 


p{uj) 




( 10 ) 


For N output photons, the photon number distribution 
is simply given by 


n{u}) = Np{uj) . (II) 

However, the extension of the wave function from one 
photon to multiple photons is somewhat conceptually 
strange. In the Heisenberg picture we discuss in the next 
section, the photon number distribution at the output is 
more naturally introduced. Also in this discussion, radi¬ 
ation pressure noise is not included. To fully take into 
account the radiation pressure noise, we need to move to 
the Heisenberg picture. 


III. WEAK VALUE AMPLIFICATION OF A 
SIGNAL IN THE HEISENBERG PICTURE 

When we deal with a radiation pressure force, it is 
convenient to work in the Heisenberg picture. To do so, 
we use the single-photon formalism for a quantum elec¬ 
tromagnetic field in an interferometer, in contrast to the 
two-photon formalism in Kimble et al. (ssj, because a 
pulsed laser has a broad spectrum around a central fre¬ 
quency. In the Heisenberg picture, an electric field is 
written as 


E{t) = E^+\t) + E^-\t) , (12) 

,13) 

£;(-) = [£:(+)]t. (14) 

Here A is the effective scattering cross-section of the 
beam, is the annihilation operator of a positive fre¬ 
quency mode w, which satisfies the following commuta¬ 
tion relations 

[at^,a,^/] = 0 , [a^i, al,] = 27rS{uj - uj') . 

Let us denote the input (incoming) and output (outgo¬ 
ing) fields at the dark port by a and b, respectively. Sup¬ 
pose that a laser pulse is injected from the bright port 
and the field is denoted by H = a + d with a classical part 
a and a quantum fluctuating part d. For brevity, we omit 
the subscript uj for the field hereafter. The input-output 
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relation for the fields is derived in Appendix and is 
given up to the linear order in vacuum fluctuations by 

& ~ { a sin cos ui^ + d sin — ia cos } , 

(15) 

where ^ = 2^ + Te and 

c\rp poo 7 

= — / —fujj{aa^ + a*a) . (16) 

m Jg 2 tt 

The definition of tq is just for convenience of notation 
and is determined so that 0/2 = ujTg{uj). The time shift 
comes from radiation pressure force exerted on a mir¬ 
ror. In the expression, m is the mirror mass and T is 
the interval of measurements. The T dependence shows 
up in merely because the mirror moves by the mo¬ 
mentum change due to the disturbance of the previous 
measurement. In the derivation, we assumed <C 1 
and neglected the higher order terms in the vacuum fields 
and w^r- 

The photon number at the output is 


of a pulse is 

^b = y ^fiuno{uj) 

= . ( 21 ) 

V Stt 

At the second line, although the range of 17 is limited to 
be 17 > —I for the probe light, we extended the integral 
range to — oo < 17 < oo because the contribution from the 
frequencies 17 < —1 is exponentially suppressed owing to 
the Gaussian tail. 

Using the definition of the phase shift due to mirror 
displacement in Eq. (|5]) and writing the phase factor in 
the signal as 

= (1 + ^)^ + 2 ’ (^^) 

we have the averaged photon-number distribution at the 
output in the form 


n = b^b 

~ lap sin^ + 2|apa;^r cos wf sinw^ 

-I- {a*d + ad^) sin^ + i{aa^ — a*a) sinw/'cosw^ . 

(17) 

Evaluating this in a coherent vacuum state |0) = |0)(i|0)Q, 
we have the average number of photons at the output 

n(uj) = (0|n|0) 

= fio(w) sin^ , (18) 


where we defined the average photon number distribution 
at the input by no{uj) = |a(a;)p. 

In what follows, we consider the Gaussian wave packet 
as is the same in the case of Schrddinger picture. How¬ 
ever, it is not a wave function of a photon but here the 
photon number distribution in the Heisenberg picture 


a(L>j) = A exp 


(w - Wo)^ 

. 


(19) 


for the frequency a; > 0. Then the spread should be 
<7i.j ^ wq. At (jj ~ 0, the distribution has an abrupt 
cutoff. This is somewhat artificial, but in most cases the 
cutoff is at the tail of the Gaussian distribution and does 
not much affect the observable signature at the output 
only if ^ wq. Defining a sideband, a; = wq + f7, and 
the normalized quantities 17 = fl/ujQ and 
the distribution is written as 


n(f7) = 


Tvr ,-0^/25'/ 




-Noe- 


sm 


2+(l + ^^)2 


(23) 


where Nq = Pq/Huio is the effective total number of pho¬ 
tons at the input. Note that the actual number of pho¬ 
ton is different from Nq since frequency is not monochro¬ 
matic. All photon number at the output is obtained by 
integrating the photon distribution over all frequencies, 
— I < 17 < oo. As well as the above integral in Eq. m, 
we extend the integral range to — oo < 17 < oo and use 
the mathematical formulas for the integrals, Eqs. m 
and (|B2D . we have the total number of photons at the 
output 


N Ki ujQ 

^ Ao 
2 


f 


27r 


n(i7) 


I — e 


' cos(0 -I- (j)) 


(24) 


The quantity o'5/)^/2 appearing in Eq. (l24ll coincides 
with the measurement strength parameter s in Eq. ®. 
When the measurement is weak (s <C I) and the post¬ 
selection is nearly orthogonal (|0| <C 1), the output pho¬ 
ton number is significantly suppressed from the input 
photon number. On the other hand, for |0| <C I, the weak 
value defined in Eq. (O is largely amplified. Therefore, 
we observe that there is a tradeoff between the magnitude 
of the signal and its statistics. 

Once we have the photon number distribution at the 
output, it is straightforward to compute the frequency 
shift and its variance. The frequency shift is 


a(I7) = A exp 



( 20 ) 


for 17 > —I and < 1. Then the average input power 


(ft) ft! Wo 
2s 


r°° dn ^ fi(f7) 

e~‘^ sin(0 -I- 4>) 

1 — cos(0 -I- (j)) 


( 25 ) 
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Figure 2. Photon number distribution at the output n{Q) in 
the unit of No/ljq. In each panel, one parameter is varied 
with fixing fiducial parameter to a^j = 1, (j) ~ 10~^ rad, and 
8 — 10“^ rad. In the upper panel, (ft = 2 x 10“^ (blue, solid), 
lO”"* (green, dotted), and 5 x 10“® rad (red, dotted-dashed). 
In the lower panel, 9 = 2 x 10“^ (blue, solid), 10“® (green, 
dotted), and 5 x lO^'^rad (red, dotted-dashed). 


changing 0 



changing 0 



Figure 3. Probability distribution at the output 
uJon(yi)/{2 tiN). The parameters for each curve in each panel 
are the same as Fig. (2] but they are normalized by the total 
photon number at the output. 


At the first line, we again approximately extended the 
integral range to from —oo. Since {Cl)(j)/2 = this 

coincides with the result in the Schrodinger picture in 
Eq. (ini). As well, the expectation value of squared fre¬ 
quency shift is 






n(0) 

N 


2s ^ 2s e ® cos{9 (ft) 
(jp' 1 — cos(6l -b (j)) 


(26) 


Again this coincides with Eq. ([71) obtained in the 
Schrodinger picture. 

Figure [2] shows the photon number distribution at the 
output n{Cl) in the unit of Nq/uiq, varying (ft and 9, re¬ 
spectively. As (ft increases, the distribution is amplified 
and more shifted. On the other hand, as 9 decreases (the 
weak value is more amplified), the distribution is more 
shifted but significantly suppressed. These behaviors in¬ 
dicate the existence of the trade-off relation between the 


amplification of the weak value and the amount of out¬ 
put statistics. In the Heisenberg formulation, the role of 
the post-selection is much clearer than the Schrodinger 
formulation. In Fig.|31 the photon number distribution is 
normalized and then is interpreted as probability distri¬ 
bution. If one see the probability distribution, the peak 
frequency is shifted without suppression of the ampli¬ 
tude. This feature often causes confusion with the am¬ 
plification of the SNR as in the Schrodinger picture. 

In this section, we considered the average photon num¬ 
ber distribution at the output, and calculated the fre¬ 
quency shift and its variance, based on the definitions in 
Eq. (l25ll and (l26l) . However, there should be vacuum fluc¬ 
tuations of an electromagnetic field, which produce fluc¬ 
tuations in the photon number and consequently quan¬ 
tum noises, i.e. shot noise and radiation pressure noise. 
To clarify the usefulness of the WVA, we need to appro¬ 
priately take into account quantum noises and evaluate 
SNR. These will be done in the next section. 
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IV. QUANTUM NOISE AND STANDARD 
QUANTUM LIMIT 


term between the vacuum fields in the same frequency 
mode but at different times by T is 


In this section, we first introduce the quantum fluc¬ 
tuations of photon number and compute shot noise and 
radiation pressure noise in a unified framework in the 
Heisenberg picture. Then defining SNR, we derive the 
SQL and the sensitivity limit to measure position of a 
mirror. 

From Eqs. dig and m, the fluctuating part of the 
photon number is 

An{uj) = niuj) — niuj) 

= {a*d + ad^) sin^ 

-I- {2|apa;^r + i{c(a^ — a*a)} sinwf cosw^ . (27) 

Then an additional contribution to the frequency shift 
due to a deviation from the average photon number is 
given by 


AQ 


Wo 


r°° dn - An(Q) 


The expectation value of AQ is zero by definition. The 
quantum noise arises from the variance of AQ, that is, 
((AQ) 2 ). 

Since in our case each frequency mode of a laser pulse 
is independent when evaluated in a vacuum state | 0 ) = 
| 0 )d| 0 )a, we have 

+ a*d}{a'{d"^)' + {a*)'d'}') = 27 rno(w)( 5 (a; — uj') , 

(28) 

((aa^ — a*a){a'{a^y — {a*)'a)') = —27rho(a;)d(w — lu ') . 

(29) 


^|aa| + a*ai| {aa^ — a*a} 

+ {aa^ — Q;*a} |aaj -|- 
= -b a*ae^‘^^} {oat - a*a} 

+ {aa^ — a*a} + a*ae*‘^^}) 

= 2i |a|^(aa^) sinwT . 

Actually there exists a nonzero correlation, which is pro¬ 
portional to sinwT. However, in practice, the frequency 
resolution of a photodetector (multichannel CCD) is fi¬ 
nite and a single frequency bin contains a number of fre¬ 
quency modes. So summing over the modes averages out 
the phase of the correlation term. This is a conspicuous 
contrast to a continuous measurement, in which the light 
spectrum is nearly monochromatic and the correlation 
term almost coherently oscillates, with the phase differ¬ 
ence of the order of ClT <C 1. Thus, we conclude that we 
can treat the fields a and oi independently in the pulsed 
measurements. 

Therefore, from Eq. dZil) and the above relations, we 
obtain 

(An(Q)An(Q')) = 27rho(a;) sin^ S{u} — w') 

+ (Cr )”o(w)rlo (lo')ujuj' sin 2uj^ sin 2u}'^ 

= 27r^d(Q - n') + w^(e^)^(Q)^(Q') , 

Wo 

(31) 

with 

/3(Q) = ho(H)(l + H) sin[(l -t- Q)()) + 0\ . (32) 


One needs to be careful about that the vacuum field a 
and that in are defined at different times, because is 
induced by a laser pulse at the past time. To distinguish 
them, we fix the subscript ” 1 ” for the vacuum field at the 
time of the previous laser pulse. Then using Eq. (HU), 42 
is given by 


The first term is shot noise, which results from the Pois- 
sonian statistics of photons and recovers Eq. (7) in (2^. 
The second term is radiation pressure noise that we orig¬ 
inally derived here in the context of the weak value am¬ 
plification in the Heisenberg picture. Then the variance 
of the frequency shift is 






— nuj[aa[ + a oi) 

ZTT 




2 rY 

m J 


dio 

27T 


[hw) no(w) 


2hu)QT 

m 




(30) 


Here we extended the integral range from — wq to —00 
owing to the exponential suppression and used Eq. (EU. 
As we prove below, it is shown that there is no correlation 
between vacuum fields at different times. The correlation 


am') = ^ j 

N 


° dn 

» 27r J 
p2~4 


dn' 

2tt 


nCl'{An{n)An{n')) 


((AD)2)sql 


(^2) sin2 e 


1 

7 + ^ 


(33) 


At the second line, we assumed a weak measurement 
( 1^1 <C 1 ) and kept the leading order term in powers 
of small (f). At the third line, we used Eq. (l30ll and the 
approximated photon number from Eq. (1241) 

iV«7Vosin2^, (34) 
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and defined 


I = Nov, 

AfvjjQGujT 
1 = - 


((Ari)^)sQL = 


COS ■ 


-2 


'\/r+^ 


2 

(jj 5 


Sin • 




I cos(0/2)|^ 


(35) 

(36) 


(37) 


These parameters, I and rj, are interpreted as the inten¬ 
sity of a measurement and the susceptibility of a mirror 
to measurement disturbance, respectively. We see in the 
variance that the first term is shot noise and the sec¬ 
ond term is radiation pressure noise and that there is 
a tradeoff between these two terms, which consequently 
gives the lower bound of total noise that one can reach. 
The minimum total noise is achievable when 7=1, for 
which the noise is ((Afl)^)sQL- This sensitivity limit is 
the so-called SQL. 

To derive minimum mirror displacement that is de¬ 
tectable by pulsed measurements, we define the SNR for 
a pointer shift: 


SNR = 


l(»)l 


Since the shift signal in Eq. (pS)) is expanded for small (j) 
as 


(f 2 ) 

the optimal SNR is 


)COt ■ 


e 


~ I 


(38) 


SNRSol = —^ \ -T 


1/2 


COS ■ 




Using 1^1 = and setting SNR = 1, we obtain the 
minimally detectable mirror displacement is 


4ql=(1 + ^) \I^ 


1/4 


'Th 


e 

COS - 
2 


- 1/2 


> 


(39) 


This exactly coincides with the standard quantum limit 
derived from an elementary argument on the SQL based 
on the Heisenberg uncertainty relation in two-time posi¬ 
tion measurements [s^. This sensitivity limit does not 
depend on the laser power of an experiment and gives the 
fundamental limit imposed by the Heisenberg inequality. 
To increase the sensitivity further, one can choose shorter 
time interval T of a measurement. However, one needs 
to increase Nq to compensate the decrease in T and keep 
7=1 {NqT = const.). Namely, to improve the sensi¬ 
tivity to £ by 10 times, one needs 10^ times larger laser 


power. This is the standard scaling appearing in quan¬ 
tum metro logy when photons in each frequency are not 
correlated [^. 

The interesting result here is that there is no ampli¬ 
fication of the SNR due to the WVA. Indeed, both the 
signal and the noise are amplified. However, in the SQL 
formula, both amplifications are canceled out and the 
weak value disappears. This can be understood intu¬ 
itively from the lower plot in Fig. [21 As 0 decreases (the 
weak value is more amplified), the distribution is sup¬ 
pressed and becomes less sharp. More strictly speaking, 
the shot noise is amplified due to the small statistics at 
the output and the radiation pressure noise is amplified 
by the same mechanism as the signal because these come 
from the displacements of a mirror. Therefore, not only 
signal but also quantum noise are amplihed. If the dis¬ 
tribution is normalized by N as in Fig. [31 one sees only 
signal amplification, which often causes the confusion in 
that the WVA is useful for improving sensitivity. 

We comment on the experimental feasibility of the SQL 
in a pulsed measurement. The SQL is given for typical 
parameters by 


£sql > 7.3 X 10-1® 


T 

10-3 s 


1/2 



(40) 


To reach the SQL, we need 7=1. For typical parameters, 
7 is written as 


7 « 1.1 X 
1 + g: 


l.Og 

m 

2 \ 1/2 
UJ 


^0 


1.8 X 1013 s-3 


e 

G^ cos - 



(41) 


Since the requirement for these parameters is not so se¬ 
vere, the order of unity for 7 and the SQL are achievable 
with state-of-art technologies. 


V. QUANTUM TECHNIQUE TO OVERCOME 
THE SQL 


There are many methods to overcome the SQL in a 
standard quantum metrology, such as using quantum en¬ 
tanglement and a squeezed state. These quantum cor¬ 
relations are also applicable to weak measurements with 
weak value amplification and enhance the sensitivity of 
parameter estimation [3^, |4^. However, what the au¬ 
thors considered was not an optical interferometry with a 
bunch of photons but just multi-time measurements and 
is not relevant to our case considered in this paper. On 
the other hand, in some cases, say, in the measurement of 
a coupling parameter between a spin-1/2 particle and co¬ 
herent light, it is shown that the Heisen berg scaling can 
be achieved with classical coherent light [dl , |42 1, though 
this does not necessarily mean that a weak measurement 
with a post-selection outperforms standard interferome¬ 
try. 




























In this section, we consider one of examples that 
circumvent the SQL in optical pulsed measurements, 
namely, the squeezed vacuum input. This is implemented 
by replacing a coherent vacuum injected from the output 
port with a squeezed vacuum. That is, the output photon 
distribution is evaluated in the state 


|0.) = |0)d®|0,)a, (42) 


where |0s)a is the squeezed vacuum for the field a at each 
frequency and is defined by 


\0s)a = S{r,,cj)s)\0)a, (43) 


with the squeezing operator 


S{rs,4>s) = exp 


duj rg^u}) 
^ 2 




(44) 


Note that the squeezing factor Vg and angle <f)g depend 
on frequency. However, for simplicity of notation, we 
omit the arguments of rg and (j)g in what follows. The 
annihilation and creation operators are converted by the 
squeezing operator as 

{rs,(j)g)aS{rg, (j)g) = a coshsinhr-s , (45) 
S\rg, (j}s)a' S{rg, (pg) = coshr^ — sinhr^, . 

(46) 


Since the squeezing operator does not change the linear¬ 
ity of a and a\ we have 

a(0,|a|0«)a = 0 a(0«|ai|0,)a = 0. (47) 


different values at different times. However, to obtain the 
analytic expression of shot noise and radiation pressure 
noise, we assumed that the squeezing factor and squeez¬ 
ing angle are constant and independent of frequency. Of 
course, in practice, the squeezing factor and angle de¬ 
pend on frequency and are not reduced to the simple 
expression in Eq. ()48|) . In that case, the quantum noise 
spectrum is somewhat degraded and becomes closer to 
the conventional spectrum without the squeezing. Since 
the vacuum fields at different times are not correlated 
due to a broadband spectrum as in the case of a coher¬ 
ent vacuum, there is again no correlation term between 
shot noise and radiation pressure noise. 

Since the frequency shift in Eq. (|38)l is the same as 
in the coherent vacuum case and does not depend on 
the squeezing parameters, it is convenient to define the 
ratio of noise in the squeezed vacuum case to that in the 
coherent vacuum case, 


^ {{AmsQL ■ 


(49) 


For specific choices of the squeezing angles, 

f±{rg,0) = ±2sinhrse=‘='''’ —1 +/± = , 

(50) 

f±{rg,Tr) = =F2sinhrse'F’'“ —>• I + /± = . 

(51) 


If 9 is small, F+ Ri /+. Then for a set of squeezing 
parameters, Vg = Vgi = rg 2 and pg = pgi = pg 2 = 0, 


R 


2 

s 


1 

2 



+ Ie 


-2ra 


(52) 


Then the quantities that is linear in a and aP, e.g. h, 
iV, An{ijj), (H), and (H^), are the same as those in the 
coherent vacuum case. However, the nonlinear terms in 
a and a^, that is, shot noise and radiation pressure noise, 
are modified. 

The derivation of quantum noise is provided in Ap¬ 
pendix [C] and the final expression when |(()| <C I is given 
by 


{{Anr)g = 


((A^)^)sql 

2 

1 + F+{,fg2, Ps2) 


+ ^ {1 + f-i'Tsl, <^sl)} 


(48) 


and for rg = = rg 2 and pg = pgi = pg 2 = tt. 


R 


2 

s 


1 

2 



+ Ie 


+ 2 rg 


(53) 


In these cases, one of the shot noise or radiation pres¬ 
sure noise are reduced and the other is enhanced. This 
means that the balance between the shot noise or radia¬ 
tion pressure noise is change and then the optimal laser 
power is also changed. However, the quantum noise never 
circumvents the SQL. If one choose the squeezing angle 
optimally as pgi = 0 and ps 2 = tt, the noise is minimized; 


R 


2 

s 



(54) 


with 

f± (rg^ps) = 2 sinh rg (sinh rg ± cos pg cosh rg) , 

0 

F+{rg,pg) = f+{rg,pg)cos'^ - . 

In the above equations, / is the same one as in Eq. (|M|) . 
We fixed the subscripts ”1” and ”2” to denote the first 
and the second laser pulses because rg and pg can take 


Therefore, we can overcome the SQL by a factor of 6“’’'* 
in SNR. In Fig. |4j we show quantum noises as a function 
of I for different choices of the squeezing angle. 

VI. CONCLUSIONS 

We have studied the WVA in position measurements in 
an optical interferometer by appropriately taking account 
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Figure 4. Ratio of quantum noise Rl with/without the 
squeezed vacuum input. The black, solid curve is conventional 
case (rjs = 0). The other curves are when the squeezing factor 
is set to rs = 1 and the squeezing angles are {(jjai, 4 >s 2 ) = (0,0) 
(green, dotted), (tt, tt) (red, dotted-dashed), and (0, tt) (blue, 
solid). 


quantum shot noise and radiation pressure noise in the 
Heisenberg picture. The Heisenberg formulation enables 
us to intuitively understand what happens when the mea¬ 
surement outcome is post-selected and the weak value is 
amplified. As discussed in Sec. im the photon number 
distribution as a function of frequency is shifted by weak 
value amplification but is significantly suppressed due to 
the post-selection. However, once this distribution is nor¬ 
malized by the total number of photons, the suppression 
of the statistics disappears. Therefore, it looks like that 
a signal is amplified without any cost and often leads to 
a wrong conclusion. However, working in the Heisenberg 
picture makes this type of confusion clear. 

Then we have took into account vacuum fluctuations 
of an electromagnetic field and derived quantum noises, 
i.e. shot noise and radiation pressure noise, in Eq. (1331) . 
By defining the SNR, we have shown that the sensitiv¬ 
ity of a mirror position measurement is limited by the 
SQL in Eq. (IMl) . which is the same as that in a stan¬ 
dard interferometry. Interestingly, the signal in Eq. (l38l) 
is amplified by the WVA but the quantum noise Eq. (1571) 
is also amplified. As a result, there is no amplification 
of SNR because in the SNR formula both amplifications 
are canceled out and the weak value disappears. We have 
also shown that the SQL can be overcome by implement¬ 
ing the squeezed vacuum input. 

Thus, we conclude that the WVA has no advantage to 
improve the sensitivity of position measurements when 
quantum noise dominates and that other quantum tech¬ 
niques must be combined with the WVA to circumvent 
the SQL. 
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Appendix A: Input - output relation 

We derive an optical input-output relation in a Michel- 
son interferometer shown in Fig. [SJ 



Figure 5. Michelson interferometer and electric fields. Here 
the initial phase offset 9 is symmetrized merely for simplicity 
of calculation. 


1. Field relation at the beam splitter 

At the beam splitter, the output field b is related with 
the other fields and Cy by 


Then using Eqs. (HU) - (HI, we have the relation 


^ 72 ■ 


(Al) 


On the other hand, the fields Cx^ Cy are expressed in terms 
of the input fields D and a as 


1 


1 


L.y — {D -|- a) , Cx — {D a) . (-^2) 
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2. Field relation in the arms 


When light goes around an interferometer arm whose 
length is L, the light takes time of a round trip in each 
arm, 2t = 2L/c, the short time delay due to small dis¬ 
placement of a mirror, 2(^ -|- 5t)/c, and a phase shifter, 
Atg for a round trip. The mirror displacement i is de¬ 
fined as differential component of mirror displacement 
between both arms because a common mode can always 
be included in the definition of L. The difference of the 
definitions between i and bi is that the former is the 
mirror displacement without radiation pressure and the 
latter is due to radiation pressure noise coming from vac¬ 
uum field fluctuations. Here Atg is defined by the phase 
factor d/2 = uj Atg{uj) and the initial phase offset 0 be¬ 
tween two arms is symmetrized merely for simplicity of 
calculation. For the other arm, the signs of £, 61, and 
Atg are reversed. 

Denoting the time delays as 

^ = 2£ + Ate , Cr = 26e , 

the electromagnetic fields departing from the beam split¬ 
ter and returning there after the reflections at the mirrors 
are connected by the relations 

- 2t + {^ + ^,)] . 

Note that is an Hermitian operator. Then using 
Eqs. dni - da, we hav the relation 

c'= , (A3) 

where the upper and lower signs correspond to y and x 
arms, respectively. 


3. Radiation pressure noise 

A pulsed light in Fig. [5] is reflected at the mirror. The 
number distributions of photons impinging on the mirror 
in each arm of an interferometer are given by 

Uy = c^Cy = -{D'^D + D + D^a) , 

Ux = cl,Cx = - D^a) . 

The number of photons included in a pulsed light is al¬ 
ways fluctuating. In the following, we ignore the higher 
order terms in vacuum fluctuations and keep up to the 
first-order terms in the amplitude of vacuum fluctuations. 
Defining D = a + dhy separating a classical part a and 
a quantum part d, we have the average photon number 


nx = n-y = |ap /2 and their fluctuations 


Any 

Aux 


- Tly fly 

= ‘^x ~ fix 

= 2 ^ced^ + a* d — aa^ 


+ a*a) 

— a* a) 


Since the classical part can be compensated exactly by 
feedback controls, we do not consider it hereafter. The 
momentum fluctuation exerted on the mirror per a pulsed 
light (twice due to reflection) is 

/■CXD 

APi = 2 -^huAmiuj) , i = x,y. (A4) 
Jo 

Let a pulse interval be T. Then at the next measure¬ 
ment time after the interval T, the mirror changes its 
position by 


= ±2(5^ 

^ 2TAPy^X 
m 

2T duj^ , . 

= ±— / — hujiaa' + a a) . (A5) 

m Jo 2 tt 

The upper sign is for the y arm and the lower sign is for 
the X arm. Also we neglect the common mode because 
it can always be included in the definition of L without 
loss of generality. 


4. Electromagnetic field at the output 

Combining Eqs. (Tni) . (IMI) . and (IA3I) . the output field 
is 


b = {iDsina;(^r + 0 + acosuj{^r + C)} ■ 

In the following, we assume |a;^r| "C 1. We also assumed 
that a is a vacuum field so that a higher order term ax 
can be neglected. Thus, plugging in D = a + d and 
neglecting the higher order terms, we obtain 

6 ~ {a sin -I- auj^r cos + d sin — ia cos } . 

(A 6 ) 


Appendix B: Integral formulas 


/ OO 

cos(26a: -|- c)dx 

-OO 

/ OO 

sin(2&a; -I- c)dx 

-OO 
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Appendix C: Derivation of quantum noise in the 
case of the squeezed vacuum input 

Since the squeezing operator does not change the lin¬ 
earity of a and o', we have 

a(0.|a|0,), =0 a(0«|a^|0,)„ = 0. (Cl) 

Then the quantities that is linear in a and a\ e.g. n, 
N, An(w), {(l), and are the same as those in the 
coherent vacuum case. However, the nonlinear terms in 
a and are modified: 

a( 0 s|a^(a')^| 0 s)a = —27r(5(a; — sinhcgcoshrs , 

a( 0 s|aa'| 0 s)a = —27r(5(a; — sinhr^coshrs , 

a( 0 s|a(a')'f| 0 s)a = 2nS{u} - w') cosh^ r* , 
a( 0 s|a'l'a'| 0 s)Q = 2TrS{oj - oj') sinh^ . 

The terms giving shot noise and radiation pressure noise 
in Eq. (EZl) are 

a( 0 ,|{aat ± a*a}{a'{ay ± (a')*a'}| 0 .)a 
= —2Trd{u} — uj') {=Fho(a;) cosh 2 rs + Re[Q;^e“*'^*] sinh 2 rs} 


As in the case of a coherent vacuum, there is no correla¬ 
tion between vacuum fields at different time because of 
a broadband spectrum. 


In the above, always appear in the product . 

Without loss of generality, we can assume that a is real 
by absorbing the phase into (j)s- So we write 


, Relate , Imfa^e 

cos (ps = -rT 5 - > sm (ps = - [—75 - 

|ar |ar 


and 


(as= 


2T 

m 


ai 
2 poo 


duj duj' 


27r 


27r 


hwhuj' 


Xai (0s|{aaJ +a*ai}{a'{a[y + (a')*a'i}|0s)c 

2 poo 




duj 

27T 


{huY 


X {|ap cosh2rsi — Re[a^e sinh2rsi} . (C2) 


In addition, we distinguish Vg and cps at the times by 
fixing the subscripts ” 1 ” and ” 2 ” for the first and the 
second laser pulses. To obtain the analytic expression of 
shot noise and radiation pressure noise, we assume that 
the squeezing factor and squeezing angle are constant and 
independent of frequency. Combining all these and using 
Eq. dUl), we obtain 


(An(w)An(a;'))s = 2TTno{oj) sin^ uj^S{uj — uj') {l -I- f+{rs 2 ,(ps 2 ) cos^ -|- (^^)sno(w)no(a;')ww'sin 2 w^sin 2 w'^ 


= - H') ^ 1 + f+{rs2,(ps2) cos^ 

UJQ 


(l + fl) 


2 2 




with 


f±{rs,(ps) = 2 sinh(sinh± cos (pg cosh Vg) 
{er)s = {er){l+f-{rslAsl)} . 


(C3) 


Then an additional term to the shot noise is 


^ dn r°^ ^ 

A2 2TT 277 


HH'27r- -6{Yi — 0')/+(rs2, (ps2) cos^ 


Wo 




dCt 


4d,,iV2 

(H 2 ) 


f+{rg2,(pg2) / ^/^"^“sin^ {l + n)(p + 6 


N 


F+{rg2,(ps2) 
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where 


- 2 

F+{rs2,<t>s2) = -^j^^Nof+{rs2,4>s2) 

X |l - - Ac^'^dl) cos[2{cj) + 0)]} 

Q 

~ /+(?'s2,</'s2)cOS^ 2 • 

At the second line, we expanded in small (f) and took the 
leading order term. The radiation pressure term is simply 
given by replacing (5^) in the calculation of a coherent 
vacuum case with 


Thus, 


{mf)s 


^ r°° ^ 


r°° dO' ~ ~ 

/ —nn'{An{n)An{n')) 

J-oo 27r 


~ {1 + F^{rs2,4>s2)} 

J\j2 ,2~4 

+ ) {1 + sin^ 0 

mmsQL 

2 


1 + F^(rs2,4 >s2) 


+ ^{1 + /-(^si, </i’si)} 


(C5) 


At the third line, we used Eqs. dMl) and dsni), expanded 
{C)s = {1 + f-irsiAsi)} iO ■ (C4) in small (j). 
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